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Jfatrix Multiplication Floating Decimal Auxiliary 

Floating Decimal Auxiliary* uses interpretive orders and 

control registers 0, 1, and 2 - 

21 without punchout — 26 with punchout 

Without punchout, 9 + 3n (2p + 2m + 11 ) + ICtop millisec. 

plus input of elements. 

With punchout, 9 + 3n (2p + 2m + 11) + 19mp milliseconds 

plus input and output time. 

Usually 9 decimals, hut depends on condition of matrices 

i»e« inversely proportional to the range of the elements in 

any one matrix. 

0* 1» 2, 3 mp words specified hy SN, n words specified "by 

SS, np words specified hy SK 

S3 - SJ 

(Address of Floating Accumulator) x 2~ 59 
(address of Code Al) x 2~™ 

Where: B is an (m hy n) matrix, A is an (n hy p) matrix, 



ACCURACY 



TEMPORARY STORAGE 



PARAMETERS 



3 
k 

5 

6 « n x 2 

7 « p x 2 

8 = mp x 2 
-39 



-39 
np x 2 ^ 

-39 



■39 



39 



e x 2 



h. is the location of the 
. c^ is the location of 



«*39 

* a °° * ! 59 

* h, x 2 yy 

-39 

* c rtrt x 2 Jy 

00 



MAXIMUM MATRICES 



DBSCRTPTICar 



and C » BA is an (m hy p) matrix* a,^ is the location 

ou 

of the first element of A* 

first element of any row i of Bi 

the first element of C * BA. e is the location of the left 

hand order to which control is transferred upon completion 

of this routine. If punch out is desired, set e = 21L 

in this routine. 

Square matrices maximum m = n=p = 2Q« Always, mp + n + 
np < 825 

This auxiliary routine performs matrix multiplication of 
two matrices B * m hy n matrix and A * n hy p matrix to 
ohtain BA * C * m hy p matrix in floating decimal form. 
The elements of C are stored in row form into successive 
locations "beginning with SH« The routine is entered hy an 
interpretive 8j order or hy an immediate 2.6 or 2k xN 
directive » 
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METHOD 



FUHCHIHG AMD SCALUffG To Insure that no element c,, of C exceeds 10 



and 



KOTB 



The Method used Is the same as that of the regular matrix 
multiplication routine, Code MO. except that no terminating 
symbol proceeds each row. 

. . — __ _ — -I* ~ - — - -o * |« Jk i « 

|h i . | should be < /10/n x 10? 2 . The elements of A and B 
are punched in row form for regular floating decimal input. 
The elements of A are punched first and then the elements 
of B. When started , the routine inputs all of A and then 
inputs B one row at a time. 

(1) If the elements in any one matrix vary over an excessive 
exponential range (> 10) the subproducts (b a ) may induce 

lO Qja. 

an error. e.g# If subproducts of an element c , are 

10* 20 , + 10 10 , + 10 15 , - 10 10 , - 10 15 , - 10?°, in that order, 

c ik will be represented by -10~ since adding first 10 

+ 10 gives 10 in the machine, and by adding lO 1 ^ - 10 10 

15 -rPO 

-10 ^ to this we get 0- Then adding 0-10 gives 

-20 
10 due to the range of the exponents exceeding 10. 

(2) If an element of A or B is represented by l/lO x 10 1 *, 
It must be remembered that l/lO has only 28 significant 
binary digits. Hence, for example, if subproducts of an 
element c^ of C are (l/lO) (4/10) x id 9 , we may not* get 
zero since the two terms in this sum will not be represented 
in the same way in the machine. The size of the error in- 
duced in such a case is proportional to the size of p* 
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3 

k 

5 
6 

7 
8 

9 
10 
11 
12 

13 
14 

15 
16 



ORDER 

kl 3F 
50 L 
26 S4 

or S5 

88 F 
OS (SK) 
03 2L 

it s8 

OK S6 
88 F 
OS (SS) 

02 4l 
2KS7 
OK S6 
05 (SS) 
&J (SK) 

84 3F 
83 3F 
8j l4L 

03 7L 

85 3F 
IS (Sfl) 
8j 16L 
12 12L 
8j S9 
22 6L 
8j 1% 
83 4l 
L5 7L 
L4 18L 
42 7L 
26 29S4 
4l 3F 
L5 7L 



Frew 3 
fcy 3 

Yesl 



i 



From 13* 
From 5 r 



1 



From 12* 
From 9 r 
By 15 



Yes 



V 




NOTES 



Clear 3 for storage 
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Enter interpretive code 

Set loop to input np elements of A» 

Input a , and store 

Is A input? 

Set loop for counting elements of C. 
Set loop to input elements of B 
Input row i of B and store. 

Is row i of B input? 

Set loop for p Z s using row i of B» 

Set loop for n products to sum. 



z V» = K(3) 

Escape to l4L 

Is E t i1 a 1k complete for column k of A? 



Store c., 
lie 

Escape to l6L 

Has C "been computed? 

Escape when finished 

Have p sums "been formed using row i of B? 

Then escape to 19L 



Increase address of a. "by p 



Heenter interpretive code 
Clear 3 for storage 



LOCATION 



HOTES 
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17 
18 

19 
20 
21 
22 

23 
2k 

25 



LO Ik 
L* k3&k 
26 15L 
00 S7 
kl S3 
L5 20L 
26' 15L 
00 SK. 
32 21L 
50 21L ' 
26 Sk 
OK S8 
05 (SBf) 
89 9F 
03 23L 
8j 25L 
OF F 
OF F 



From 13 



From 12 



From 2k 
By 2k 
Yes 

From 2k* 



Reset address of a., for next sumnatiqa 



Address used 

Clear S3 for 83 order. 



Address used. 
For Print out 

Enter Code Al 

Set cycle for np elements of C» 

c ik 

Print c., to 9 decimals 
lk 

Bare mp elements of C been printed? 
Escape to stop. 



